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^ : ABSTRACT 

rr We examine effects on the charge states of energetic ions associated with gradual solar 

flares due to shock heating and stripping at high ion velocities. Recent measurements 
of the mean charges of various elements after the flares of 1992 Oct 30 and 1992 Nov 2 
allow one to place limits on the product of the electron density times the acceleration or 
coronal residence time. In particular, any residence in coronal loops must be for <0.03 
s, which rules out models of coronal transport in loops, such as the bird cage model. 
The results do not contradict models of shock acceleration of energetic ions from coronal 

(N| ■ plasma at various solar longitudes. 
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Introduction 



particles are freshly accelerated on open field lines at 



Although particles from solar flares have been ob- 



served for 50 years ( Forbush f946 ), there are still 



many unresolved questions about their acceleration 
sites and mechanisms. Much recent research has fo- 
cused on a simple classification scheme proposed by 
Pallavicini, Serio, & Vaiana (1977). In modern termi- 
nology, "impulsive" solar flares are typically defined 
as those with a short (<; 1 h) duration of X-ray emis- 
sion, while "gradual" flares have a longer X-ray dura- 
tion. These two classes of flares have been found to 
have several distinguishing characteristics. For exam- 
ple, X-ray emission from impulsive flares is observed 
to come from compact regions at low coronal heights 
(<; 10 4 km), while X-rays from gradual flares tend to 
arise from broader regions or from large coronal loops 
up to ~ 10 5 km above the photosphere. In addition, 
it is widely believed that the main acceleration mech- 
anism for particles escaping from impulsive flares is 
stochastic acceleration (second-order Fermi acceler- 
ation by gyroresonant plasma waves; e.g., Temcrin 



Sz Roth 1994 pVIiller fc Vinas 1993J ), while gradual 
flares are associated with shock acceleration higher in 



the corona (e.g., Cane, McGuire, & von Roscnvingc 



1986, Lee fc Ryan 1986 ]), th ough not all observ ations 
support this paradigm (e.g., Mazur et al. 1992 ). 

A more controversial issue is why solar cosmic rays 
are observed at locations that are magnetically con- 
nected to solar longitudes far from the longitude of the 
flare site. Because interplanetary diffusion is highly 
anisotropic, inhibiting motion perpendicular to the 
magnetic field (e.g., Palmer 1982 ), for decades it was 
assumed that particles are transported to other so- 
lar longitudes within the solar corona, whence they 



escape to travel along the interplanetary magnetic 



field to the observer. The seminal quantitative model 
of Rcid (1964) assumes isotropic, two-dimensional, 
coronal diffusion, and gives specific predictions for 
the injection rate of particles into the interplanetary 
medium as a function of time and coronal distance 
from the flare site. It is not clear exactly what mech- 
anisms give rise to coronal diffusion, though Newkirk 
& Wentzel (1978) presented the "bird cage" model 
of rigidity- independent coronal transport, in which 
flare-accelerated particles bounce back and forth in- 
side coronal loops and occasionally transfer to other 
loops at the footpoints. On the other hand, it has 
recently been proposed that for gradual flares, with 
are often associated with large, interplanetary shocks, 



different hcliolongitudes ( Mason, Glocckler, fc Hoves- 
tadt 1984, Rcamcs 1990| ). It is further assumed that 



coronal transport does not occur, and that for im- 
pulsive flares, which are usually not associated with 
interplanetary shocks, the narrower longitudinal dis- 



persion (Reames, Cane, & von Rosenvinge 1990) is 
due to the spreading of magnetic field lines in the so- 
lar corona or the interplanetary medium. Therefore, 
it is no longer clear whether several decades worth 
of solar cosmic ray observations were providing in- 
formation about the solar flare site, as previously as- 
sumed, or about acceleration at coronal or interplane- 
tary shocks. Further progress in the interpretation of 
the information embodied in solar cosmic ray obser- 
vations thus depends on the resolution of this issue. 

Recent observations which can shed much light on 
this issue concern the charge state distributions of 
ions accelerated as a result of solar flares. Such dis- 
tributions for various elements provide a rich source 
of information on the conditions of particle acceler- 
ation and escape from the corona. Results of the 
ULEZEQ instrument on board the ISEE-3 spacecraft 
indicated that for several gradual flares, the charge 
states of 9 different elements were not consistent with 



the same temperature (Luhn et al. 1984), but rather 
had apparent ionization temperatures of 1 to 8xl0 6 
K. In contrast, summed results for several less pow- 
erful impulsive flares suggested that charge states of 
Si and Fe were characteristic of significantly hotter 
plasma, as might be expected if particle acceleration 
takes place at the site of a compact, impulsive flare 
(Luhn et al. 1987). More recently, measurements 



by three instruments on board the SAMPEX mis- 



sion (Octlikcr et al. 1995, Leske et al. 1995, Mason et 



al. 1995) have generally confirmed the earlier results 



for gradual flares, and also provide measurements for 
more elements, with better statistical accuracy, and 
for a broader energy range. A possible explanation for 
why the mean charges of energetic ions from gradual 
events are not consistent with a single temperature 
was given by Mullan and Waldron (1986), who pro- 
posed that photoionization by flare X-rays changes 
the ionization equilibrium in the plasma from which 
ions are accelerated, and obtained good quantitative 
agreement with the apparent ionization temperatures 
for various elements. Another possible explanation 
in terms of freeze-out temperatures was discussed by 
Mason et al. (1995). 

In this report, we examine effects on the charge 



states of energetic ions associated with gradual solar 
flares due two processes which would be expected to 
occur at the acceleration site. Although the physical 
conditions of the ion acceleration are poorly known, it 
is possible to roughly estimate these effects. If we ac- 
cept that ions escaping after gradual flares arise from 
shock acceleration high in the corona, then one might 
expect significant electron heating, which is a gen- 
eral phenomenon of shocks in astrophysical plasmas 
(Schwarz et al. 1988). Another effect is that as the 



ions' velocity increases, their electrons arc more likely 



to he stripped ott ( [riovestadt ct al. 1984| , |Lcskc ct al 



19E|3[^ Wc point out tnat it tnc ions spend a sul- 
ficiently long time at the acceleration site or in the 
solar corona, then these effects would influence the 
charge states in a manner that is inconsistent with 
the observations. This places limits on the accelera- 
tion time, and on the time spent in the solar corona 
after acceleration. These limits in turn rule out mod- 
els of ion acceleration or coronal transport (such as 
the bird cage model) which require that escaping ions 
spend over ~ 0.03 s in a magnetic loop, while they do 
not contradict models in which ions promptly escape 
after acceleration on open field lines. 

2. Effect of Shock Heating and Ion Velocities 
on Mean Charge States 

As mentioned above, there are several recent re- 
ports of high-quality data on the mean charges of so- 
lar cosmic ray ions. Here we will mainly consider 
the results of Leske et al. (1995), who measured the 
mean charges of 12 elements with a high statistical 
accuracy at the relatively high energy of ^15-70 McV 
nuclcon -1 . Very similar results were obtained for two 
large, gradual flares on 1992 Oct 30 (X-ray magnitude 
XI, importance 2B, at 22° S 61° W) and 1992 Nov 
2JX9, 2B, at 23°S, 90°W; |Solar-Geophysical Data 
1991) , 



to place conservative limits on the acceleration or 



coronal residence time for solar cosmic rays, we as- 
sume that the charge states are originally characteris- 
tic of a relatively cool coronal plasma in thermal equi- 
librium, and that other effects, such as shock heating 
or increased stripping for faster ions, would tend to 
increase the charge. The lowest apparent ionization 
temperature that was accurately measured by Leske 
et al. (1995) is that for Si, 1.75±0.11xl0 6 K, which 
was essentially the same on both 1992 Oct 30 and 
1992 Nov 2. For a conservative limit, we use an up- 



stream plasma temperature of T u = 1.5xl0 6 K, which 
is also the average value obtained by Mullan & Wal- 
dron (1986) for their explanation in terms of X-ray 
photoionization. 

The first effect which we consider is that of electron 
heating by shocks propagating through the corona, 
since it is widely believed that cosmic rays asso- 
ciated with gradual flares are accelerated by such 
shocks. Observations of a large number of fast- 
mode shock crossings of interplanetary shocks and 
bow shocks of Earth, Jupiter, and Uranus (Thomson 
et al. 1987; Schwarz et al. 1988|) imply that a fixed 



fraction of 0.12±0.05 of the incident ion ram kinetic 
energy goes into electron heating, i.e., into increasing 
the electron enthalpy. From Figure 3 of Schwarz et 
al. (1988), we see that the jump in the electron tem- 
perature, AT e , across the shocks in their sample is 
typically given by 



AT e = (3 to 9K) x 



V 2 - V? 

v u r d 



km 



2 -2 
S z 



(1) 



Vf > 2.5 x 10 5 km 2 



especially for V 2 
V u and Vd are the bulk flow speeds upstream and 
downstream, respectively, in the shock (de Hoffmann- 
Teller) frame. 

To estimate the shock speed, we note that mag- 
netohydrodynamic simulations indicate that the ini- 
tial coronal mass ejection (CME) speed is roughly 
twice the average propagation speed of the forward 
shock to the Earth (Hcras et al. 1995). Based on 



the K p index of geomagnetic activity (U.S. National 
Geophysical Data Center, via World Wide Web), we 
estimate the initial CME speed for both flare events 
to be V u ~ 1500 km s _1 . Assuming Vd — U u /4 for a 
strong shock, the downstream electron temperature, 
Td, is estimated to be 6 to 19xl0 6 K. Note that since 
shock speeds determined from Type II radio bursts 
can be high er than initial CME speeds (e.g., Heras 
ct al. 1995 ), using such speeds would yield higher 



downstream temperatures. A conservative limit on 
nt is obtained by considering Td = 7 x 10 6 K, which 
was also invoked in a mechanism proposed by Luhn 
et al. (1984) in which ions are suddenly heated to 
this temperature and the different apparent ioniza- 
tion temperatures for various elements are ascribed 
to different rates of approach to the new ionization 
equilibrium. For comparison, we also consider the ef- 
fect of a higher value, Td — 1.5 x 10 7 K. 

To calculate charge state distributions, we used the 
ionization and recombination rate coefficients of Shull 



& Van Steenberg (1982), correcting the misprints in 
those tables reported in the errata (ApJS, 49, 351) 
and by Arnaud & Rothenflug (1985). The rate of 
change of n q , the number density of ions of charge q, 
is given by 

d/n 

— - = n e [I q - 1 n q - 1 - (I q + R q -i)n q + R q n q+1 ], (2) 

where I q and R q are temperature-dependent rate co- 
efficients (in cm 3 s _1 ) of ionization from charge q to 
q + 1 and recombination from q + 1 to q, respectively. 
For ten elements, and different nt values, we used Eu- 
ler's method to calculate the mean charges expected 
due to a) shock acceleration, in which particles are 
sometimes at T u and sometimes at T^ (we used equal 
times for a conservative limit), and b) the mechanism 
of Luhn et al. (1984), in which ions are immersed com- 
pletely in a new Tj. Table 1 shows the mean charges 
of various elements for case a) . Note that for a higher 
2d, smaller values of nt yield similar results. For case 
b), almost identical mean charges were obtained for 
nt values half as large. It is clear that neither mecha- 
nism is able to explain the observations. In particular, 
at nt = 2 x 10 10 cm~ 3 s, the mean charges calculated 
for Si, S, Ar, Ca, and Ni were all above the la limits 
of Leske et al. (1995), while those for O, Ne, and Fe 
were too low. We therefore place a conservative limit 
of nt < 2 x 10 10 based on shock heating, where t is 
the total residence time at the shock. 

Next we consider the effect of the ion velocity. The 
new observations of solar cosmic-ray charge states af- 
ter gradual flares at higher energie s, ;> 10 MeV ( Leske 
et |il. 1995| , petliker et al. 1995] , |Tylka et al. 1995[ ), 
are in general agreement with those at lower energies 
(Luhn ct al. 1984, Mason et al. 1995; note, however, 
that different mean charges are obtained for Fe, as 
discussed in the latter reference), which implies that 
little or no additional stripping occurs during the ac- 
celeration to higher energies or escape from the corona 
( Leske ct al. 199q ). Firstly, let us consider a limit on 
the integrated electron density times residence time in 
the corona after ions have been accelerated to their fi- 
nal energy, assuming that all changes in charge states 
relative to equilibrium at the plasma temperature 
occur during post-acceleration motion through the 
corona. Ions of 15 MeV nucleon -1 (at the lower end of 



the energy range measured by Leske et al. 1995 ) move 
at « 6.4 times the rms thermal speed of electrons at 
our assumed plasma temperature of T u — 1.5 x 10 6 
K. At these relative speeds, collisional ionization is 



far above threshold, and varies slowly with the rela- 
tive speed, so we can neglect the effect of the thermal 
spread and approximate the ionization cross sections 
by assuming w rcl = v ion . We use equation (2) to cal- 
culate the charge state distributions, but neglect the 
recombination rates because the charge states char- 
acteristic of T u are far from the equilibrium for such 
a high velocity. For ionization rates, we employ the 
energy-dependent ionization cross sections a q of Ar- 
naud and Rothenflug (1985), and set I q = a q v le \. 
Again, the charge state observations for various el- 
ements cannot be explained for any value of nt. As 
for the case of shock heating, the tightest limits on 
nt are obtained for Ar, Ca, S, and Si. The nt val- 
ues for which the calculated mean charges exceed the 
observed mean charges by one standard deviation are 
shown in the third column of Table 2. From these 
we conclude that nt should not exceed about 3xl0 9 
cm -3 s during transport through the corona. 

Similarly, if we assume that all deviations from 
thermal-equilibrium charge states occur during ac- 
celeration, we can derive limits on the acceleration 
time. Here we need to consider the changing veloc- 
ity of the ion. In fact, for slower ions the rate of 
ionization is faster, though the ultimate equilibrium 
charges would be lower. To be conservative, we only 
consider the final e-folding in momentum from 61 to 
167 MeV c , corresponding to kinetic energies from 
2 to 15 MeV, and assume as a rough approximation 
that the rate of change of log(p) is constant (which 
requires that the spatial diffusion coefficient should 
be roughly constant over that energy range) . The re- 
sulting nt values for which calculated mean charges 
exceed the observed mean charges by one standard 
deviation are shown in the second column of Table 
2, where t represents the momentum e-folding time. 
We can conclude that the maximum momentum e- 
folding time during acceleration, which is commonly 
called the acceleration time, is given by nt < 2 x 10 9 
cm~ 3 s. Since the limits we derive for acceleration and 
coronal residence times are of a similar magnitude, we 
can alternatively say that nt < 3 x 10 9 cm" 3 , where 
t is the acceleration time plus the residence time in 
the corona. The earlier limit derived from the effect 
of shock heating is less stringent than that from the 
effect of the ion velocity, so we will use the latter in 
the following discussion. 



3. Discussion 

We have implicitly made the standard assumption 
that the charge states of solar cosmic rays are fixed 
in the corona and unaffected by interplanetary trans- 
port. This is justified by considering the above lim- 
its. Given tha t solar wind densities are on the order 
of 4-12 cm~ 3 (Foukal 1990) near Earth, and vary as 



the inverse square of the radius from the Sun, ions 
could travel for years throughout the inner solar sys- 
tem without being substantially stripped. However, 
for the much higher densities within the solar corona, 
these limits on nt correspond to much shorter times. 

A related point is that Tylka et al. (1995) have 
stressed that most measurements of the mean charge 
of iron are incompatible with acceleration directly 
from the ambient solar wind (Reames 1990). This 



poiiit is supported by previous theoretical and obser- 



vational studies (e.g., icrasawa fc Scholcr l!J8!j| 7 Tan 



ct ; d. 1990] ), which found that at energies ^ 1 McV, 
interplanetary shocks only re-accelerate an existing 
population of somewhat lower energy. The observed 
charge states seem to require an initial acceleration 
from coronal material. This does not rule out the 
possibility that a CME-driven shock accelerates par- 
ticles while inside the corona, or perhaps near the 
Sun (Kahlcr, Reames, fc Sheeley 1990|) while pusl 



ing coronal m aterial in front of it (Boberg, Tylka, & 
Ad jims 1996J) . 

We consider the implications of limits on the prod- 



uct ofTHiFcnxictronTIcli^^ 



or residence in the corona for models of coronal trans- 
port. While the flare of 1992 Oct 30 was at a solar 
longitude of 61° W that could be magnetically well- 
connected to the detectors near the Earth, the flare 
of 1992 Nov 2 was at « 90° W, so the hypothesis 
of coronal transport would have to account for a net 
displacement of about 4 x 10 5 km. If we consider the 
limit of nt < 3 x 10 9 cm -3 s in the corona after and 
during acceleration, the electron densities in various 
regions of the corona (and the corresponding limits on 
t) are roughly 10 11 " 12 cm" 3 (0.003-0.03 s) in coronal 
loops, and 5 x 10 7 cm~ 3 (60 s) at a height of 10 5 km 
in the quiet Sun (Foukal 1990). As field lines outside 
the loops rapidly open up with height, the density 
rapidly decreases, until n ~ 10 5 cm~ 3 (t^3x 10 4 s) 
at r = 2Rq. 

These limits essentially rule out the "bird cage" 
model of rigidity- independent coronal transport (New- 
kirk & Wentzel 1978) for escaping ions from these 



gradual flares. In this model, particles are stored in- 
side coronal loops, mirroring back and forth between 
the footpoints and occasionally transferring to a dif- 
ferent loop. However, the above limits show that 
the ions which escape could have traveled at most 
1.6xl0 3 km inside a coronal loop, which would only 
be a small fraction of the loop's length, and is far 
shorter than the coronal distance from the site of the 
1992 Nov 2 flare to the footpoint of the detectors' 
magnetic field lines. Thus if the ions which are ob- 
served are accelerated inside a loop, the acceleration 
time must be very fast (< 0.02 s), and the ions must 
escape immediately and cannot be transported along 
the loop. 

If a mechanism for coronal transport does not 
involve the motion of ions in cor onal loops, e.g., 
by tr ansport along current sheets ( Fisk fc Schattcn 
1972| ; note that such transport would be rigidity- 
dependent), the above limit on the time spent at 
a coronal height of ~ 10 5 km allows a 15 McV 
nucleon -1 ion to travel about 3xl0 6 km. However, 
the above limits on nt are very conservative, in that 
they permit an ion's mean charge to change from the 
value characteristic of 1.5 x 10 6 K to the measured 
values, while in fact there seems to be little change in 
the mean charges of observations spanning two orders 
of magnitude in energy. It is possible that SAMPEX 
observations could yield d(Q)/dE for a single instru- 
ment, placing more severe constraints on nt and on 



coronal transport models in general (see also Octliker 



ct al. 1997) 



Our results favor models in which the acceleration 
of ions observed after gradual flares takes place on 
open magnetic field lines (e.g., Lee fc Ryan 1986] ), 
and in which the longitudinal distribution of solar 
cosmic rays from such events is due to acceleration 
distributed over a range of heliolongitudes. 

Since the effects of shock heating and ion velocities 
cannot explain the different apparent ionization tem- 
peratures for various elements, the best quantitative 
explanation is that of Mullan and Waldron (1986). 
An requirement of this model is that similar X-ray 
fluxes must be present at the acceleration site for all 
of the observed flares, over a time sufficient for ion- 
ization equilibrium to be established, before a shock 
arrives and accelerates the ions. Note that the X-ray 
flux should decrease with increasing coronal distance, 
and be eclipsed for low heights at large coronal dis- 
tances. If one believes that observed ions are accel- 
erated near the footpoint of the observer's field line, 



then the fact that there is no apparent dependence of 
the mean charges on the coronal distance (e.g., for the 
flares of 1992 Oct 30 and 1992 Nov 2) may be prob- 
lematic for either the distributed acceleration model 
or the photoionization model. 
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TABLE 1 
Mean Charges due to Shock Heating 3 





Element Measured 13 


T c i — T u 


T d = 7 x 10 6 K 


T d = 1.5 


x 10 7 K 




nt = 1 x 10 10 nt = 2x 10 10 


nt = 5 x 10 9 


nt = lx 10 10 


l> 


N 6.47±0.20 


5.82 


6.09 6.28 


6.12 


6.33 


On 
On 


6.95±0.20 


6.20 


6.42 6.60 


6.48 


6.71 


Ne 8.53±0.27 


7.99 


8.05 8.10 


8.10 


8.19 




Mg 10.30±0.34 


9.89 


9.96 9.99 


9.98 


10.03 


Si 10.54±0.37 


9.65 


10.65 11.14 


10.64 


11.15 


P-H 


S 10.84±0.44 


9.19 


10.71 11.42 


10.76 


11.57 


OS 


Ar 10.08±0.91 


9.54 


11.02 11.81 


11.13 


12.05 


Ca 11.46±0.49 


10.35 


11.48 12.18 


11.66 


12.52 




Fe 15.18±0.73 


10.39 


12.68 13.09 


12.64 


13.19 


> 

o 


Ni 12.62±1.30 


9.32 


12.70 13.76 


12.52 


13.92 












Note. — Units of nt are cm 
a Calculated assuming that 


" 3 s. 
ions spend equal times at T u — 1.5 x 10 6 


K and at T d . Italicized values are 


O 


those 

in excess of the measured values plus 
b Leske et al. 1995. 


one standard deviation. 
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TABLE 2 










Limits on nt for Fast Ions 







Element 



nt during acceleration, 

2 to 15 MeV nucleoli" la 

(10 9 cm" 3 s) 



nt after acceleration, 
15 MeV nucleon -1 



(10 9 



cm 



Si 

S 

Ar 

Ca 



2.6 
2.3 
1.2 
1.6 



4.8 
4.4 
2.2 

2.6 



a Over a factor of e in momentum, i.e., from 61 to 167 MeV 



nucleon 



